Towards a quantum theory of chiral magnetic effect 
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We discuss three possible ways of addressing quantum physics behind chiral magnetic effect and 
electric charge fluctuation patterns in heavy ion collisions. The first one makes use of P-parity 
violation probed by local order parameters, the second considers CME in quantum measurement 
theory framework and the third way is to study P-odd * P-odd contributions to P-even observables. 
In the latter approach relevant form-factor is constructed and computed for weak magnetic field in 
confinement region and for free quarks in strong field. It is shown that the effect is negligible in 
the former case. We also discuss saturation effect - charge fluctuation asymmetry for free fermions 
reaches constant value at asymptotically large fields. 

PACS numbers: 12.38.Aw, 12.38.Mh 



I. INTRODUCTION 

One of the main theoretical challenges of modern 
quantum chromodynamics (QCD) is to build a detailed 
theoretical picture of strong interaction physics relevant 
for heavy ion collisions. Currently running experimental 
programs have already brought lots of exciting results. 
Despite tremendous progress in understanding, rich pat- 
tern of observed effects is still waiting for being placed 
into coherent theoretical picture based on QCD. 

In the course of studies of hadronic matter at large 
temperatures and/or densities one can make use of the 
scale separation allowing to neglect effects of weak and 
electromagnetic interactions in most cases. A possible 
interesting exception is pointed out in [H, . When rel- 
ativistic ions undergo noncentral collision, strong mag- 
netic field is generated in the collision region. The 
typical magnitude of this field is estimated as \J eB = 
10 100 MeV, i.e. of the order of dynamical QCD 
scale. Correspondingly, any studies of strongly inter- 
acting matter in heavy ion collisions have to take the 
effects of this abelian magnetic field into account. Of 
particular interest in this respect is the so called Chiral 
Magnetic Effect (CME). The physics behind it can be 
explained in several different but complementary ways 
[H-flEj]. Let us consider nonzero density of one flavor 
of free massless quarks in external magnetic field. Sup- 
pose there are unequal chemical potentials for left and 
right handed quarks: (J.l 7^ fJ-R- When it can be shown 
that a nonzero classical electric current flows along the 
magnetic field (see [|[ and references therein, see also 
[26[ for another prospective): 
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where 2/15 = fJ>R — [il- The physical reason for this chi- 
ral charge excess to electric charge current conversion is 
quark magnetic moment interaction with the magnetic 
field (which is of different sign for positively and nega- 
tively charged quarks) together with the correlation of 
spin and momentum for chiral fermions. Both sides of 



([T]) have of course the same transformation properties 
under P- and CP-parity conjugation. Many different 
aspects of CME have been extensively discussed in the 
literature and there is no doubt that CME is a robust 
theoretical effect. However it is not a simple task to ap- 
ply this clear physical picture to real processes described 
by nonperturbative QCD. One of the most important 
questions on this way is about physical origin of chiral 
chemical potential /15, which is absent in fundamental 
QCD Lagrangian. In original picture [7[ appearance 
of effective ^5 ^ is a nonperturbative QCD effect, 
caused by interaction of quarks with topologically non- 
trivial gluon field configurations above the phase tran- 
sition. The physical explanation goes as follows. As 
is well known the topological charge in the QCD vac- 
uum fluctuates as described by Veneziano-Witten for- 
mula m m 



(2) 



X = / d 4 x(GG{x)GG(0)) cx F l m 1 r)l 



where the nonperturbative parameter in the r.h.s. scales 
as A.q CD which means that topological charge fluctuates 
over Euclidean 4- volumes of typical size determined by 
nonperturbative QCD scale. It is worth stressing that 
these fluctuations are quantum, i.e. the states of differ- 
ent topological charge are to be summed over for what- 
ever Euclidean 4- volume V and one always has 



d i x(GG(x)) = 



(3) 



In other words, ^ vanishes because local average 
(GG(x)) — and not due to the presence of integra- 
tion over the volume V. There is no special space-time 
fluctuation pattern in the problem other than the cor- 
relator @ (and higher ones). 

The situation however may change at nonzero tem- 
perature/density. Since the Euclidean 0(4) invariance 
of the vacuum is broken in this case, one can think of 
different fluctuation patterns in spatial and in temporal 
directions. Moreover, since in real collision experiments 
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external conditions are time-dependent they can play a 
dual role of the background and of a measuring device. 
In other words the meaning of averaging in ((3]) changes: 
one has to integrate only over those field excitations 
which are present at a given Minkowski 3-volume for a 
given time period and the problem becomes essentially 
non-stationary in this sense. One can say that the aver- 
age over fields (..) becomes 1^-dependent. Such quantity 

- physically corresponding to a "single event" - can in 
principle be non- vanishing. Of course it is natural to 
expect that random character of fluctuations leads to 
zero result for ([3]) after averaging over many events. 

The CME is often considered as a reasonable expla- 
nation of outgoing particles electric charge asymm etry 
observed at Relativistic Heavy Ion Collider (RHIC) [29| 

- in ^/sWn = 200 GeV Au+Au and Cu+Cu colli- 
sions. The latter effect can be described as follows. For 
noncentral collision one can fix the reaction plane by two 
vectors: beam momentum and impact parameter (with- 
out loss of generality this is always chosen as 12 plane in 
the present paper and no adjustment angle HfRp is in- 
troduced) . Thus angular momentum of the beams (and 
the corresponding magnetic field) is oriented along the 
axis 3. The azimuthal angle cj> £ [0, 2tt) is defined in the 
plane 23. With this notation, in any particular event 
one studies charged particles distribution in <j> using the 
following conventional parametrization 

— oc l+2v\ ± cos0+2«2 ± cos2</>+2a± sin </>+... (4) 

dcj) 

The coefficients Vi t ± and «2,± account for the so called 
directed and elliptic flow. They are believed to be uni- 
versal for positively and negatively charged particles 
with good accuracy. The coefficients a + and a_ de- 
scribe charge flow along the third axis, i.e. normal to 
the reaction plane. This P-parity forbidden correlation 
between a polar vector (electric current) and the axial 
one (angular momentum) is considered as a signature of 
P-parity violation in a given event with a± ^= 0. On the 
other hand, the random nature of the process dictates 
(a+) e — ( a -)e = (there the averaging over events is 
taken) . 

Trying to construct a theory of the phenomenon one 
has first to choose adequate language. Since at the end 
the heavy-ion collision is a scattering problem, the ulti- 
mate framework would be S-matrix and inelastic scat- 
tering amplitudes formalism with two colliding ions as 
incoming particles. Due to extreme complexity this way 
seems to be totally hopeless. Instead one uses some ef- 
fective theories like hydrodynamics to predict distribu- 
tion of outgoing particles. In the particular problem 
of charge fluctuations asymmetry the crucial point dis- 
tinguishing different theoretical models is whether the 
currents of interest are treated as classical or as quan- 
tum. In the former case one makes use of the expres- 
sion ([T]) as classical equation. The quantum nature of 
the problem here is hidden in a theory for and cor- 
responding correlators and fluctuations for this effec- 



tive chiral chemical potential. In the later case one is 
to consider quantum averages like (f2|j M |f2), (f2|j M j„|f2) 
etc. and to understand (H} as operator relation. How- 
ever if one takes diagonal matrix element of (|TJ) in the 
vacuum the answer is of course trivial: (0|j|0) = even 
for nonzero external magnetic field. The absence of net 
electric current is directly related to the fact that fun- 
damental QCD Lagrangian contains no such quantities 
as fi L or fi R . 

We discuss three basic complementary ways to ad- 
dress quantum nature of CME in this paper: 

1. To make use of P-parity violation probed by local 
order parameters 

2. To consider CME in quantum measurement the- 
ory framework 

3. To study P-odd x P-odd contributions to P-even 
observables. 

We discuss all these approaches in the present paper 
and start with the first one in the next Section which is 
phenomenologically the simplest. 



II. P-PARITY VIOLATION PROBED BY 
LOCAL ORDER PARAMETERS 

As is well known quantum field theoretical averages 
of local operators have typically the following leading 
contribution: 

(n\o{x)\n) occ- A do (5) 

where A is ultraviolet cutoff and numerical constant c is 
generally non- vanishing if c = is not protected by some 
symmetry. Therefore the crucial step in the discussed 
problem is to model transition from local microscopic 
current j M to nonlocal macroscopic one J^. It is done 
by taking the matrix elements of the current over the 
medium degrees of freedom |$) from full state vector 
|n> = |*)®M: 

3n(x) = <H^O) 'h « ( $ l J dx Pv{x)j^(x)\<^>) (6) 

Here the function py (x) defines the measure of integra- 
tion over "physically infinitesimal volume", as is usual 
in condensed matter physics. 

The second important ingredient is the existence of 
the medium itself. For phenomenological purposes it is 
not important what particular kind of microscopic de- 
scription for the medium is chosen. What does matter 
is Lorentz symmetry breaking following from the exis- 
tence of a distinguished frame - the medium rest frame. 
In the simplest cases of uniform medium characterized 
by nonzero temperature/density it is usually parameter- 
ized by a unit vector - the medium four-velocity, so 
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that for applied uniform electromagnetic field one has 
the standard text-book answer for induced current 

(4>\ J^) ex u v F^ (7) 

We say about local parity violation in a state |f2) 
when a local parity-odd operator O(x) = —PO(x)P^ 
has nonzero expectation value in this state 

(fi|0(a?)|n) + (8) 

for example (ip^j^ip) ^ 0. The condition of locality here 
is important. Operationally it means that the opera- 
tors and their products are defined at the scale a ~ A -1 
where A is ultraviolet cutoff. For nonlocal averages, on 
the other hand, it is not a problem to have nonzero 
P-odd matrix element, e.g. (jo(x)j3(y)}- The medium, 
characterized by finite coherence length, brings physical 
meaning to this nonlocality. For example, in a medium 
with applied uniform electromagnetic field nothing for- 
bids to have P-odd correlation between axial current 
divergence and the vector current: 

(0|j M aj 5 l^) km"^ (9) 

where F — -f aF a P 

To feel the physical meaning of @ let us imagine 
radial distribution of velocities v of the matter in a uni- 
form magnetic field B. If the divergence d J 5 is also 
uniform in the ("fireball") volume, the charge density 
is to be of different sign above and below the reaction 
plane: 

((j>\J dJ 5 \<p) cx v-B (10) 

In medium rest frame characterized by — (1,0,0,0) 
for uniform magnetic background, the electric current J 
flows along the magnetic field B. 

It seems quite natural to interpret (JSJ in the following 
way: as soon as the concept of a medium can be ap- 
plied to the discussed problem one can easily construct 
classical nonzero local P-odd parameters without spec- 
ifying any particular " chiral microscopy" . The medium 
(manifested by existence of the selected frame) is cru- 
cial in two aspects: first, it allows to consider meaning- 
ful local objects and not badly divergent quantities like 
([5]) and second, by Lorentz invariance breaking it pro- 
vides invariant meaning for the electric and magnetic 
fields, thus making possible correlations between local 
(in macroscopic sense!) operators of different parities. 
We also see here the importance of the uniformity con- 
dition: if d J 5 is short-correlated, there is no net effect. 
This brings us back to the question about dynamical 
scales hierarchy. 

There is a deep question behind the above consider- 
ation: if the (microscopic) current non-conservation is 
anomalous (e.g. in ©) - how is this fact encoded in 
equations for macroscopic, effective currents? We leave 
aside the discussion of this important point and refer an 



interested reader to [4l| where this question is addressed 
in hydrodynamic setup. 

From heavy ions collision point of view the P-parity 
violating average © is not an observable by itself. The 
reason is physically clear: instead of measuring com- 
ponents of vector (electric) current and axial (chiral) 
current and studying their correlation, only the quanti- 
ties of the former kind are being measured (in the form 
of final particles electric charge distribution) . As for the 
latter quantities related to chiral charge - it is assumed 
that the quark-gluon medium created after collision of 
two heavy ions plays itself a role of the measuring device 
performing effective measurement of topological charge 
in the corresponding space-time region. It should be 
mentioned that this is a rather strong assumption: to 
say that one part of some quantum system can mea- 
sure (in classical sense) the state of another part of the 
same system means in fact to address some scenario 
for decoherence of the subsystems and information loss. 
To model this effect one has to adopt the language of 
quantum measurement theory. This is done in the next 
Section. 



III. CME IN QUANTUM MEASUREMENT 
THEORY FRAMEWORK 

It is possible to understand (JTJ) as a correlation be- 
tween preferred direction of outgoing electric charge dis- 
tribution asymmetry and the magnetic field in a partic- 
ular event. The sign of this P-parity odd asymmetry 
is fixed by the sign of effective ji^ in this event (and 
of course varies randomly from event to event due to 
topological neutrality of QCD vacuum). The quantita- 
tive theory would require information about distribu- 
tion function of effective 

A simple quantum-mechanical analogy can be useful 
to illustrate the point. In one-dimensional bound state 
problem with P-parity even potential V(x) = V(—x) 
one has (x) — J xdx\tpo(x, t)\ 2 = where ipo(x,t) is 
the ground state P-parity even wave function. On the 
other hand, performing particle position measurements 
on ensemble of N identical systems all in the ground 
state one gets sequence of positive and negative num- 
bers X\, X2, ■■, %n (with some uncertainties determined 
by the measuring device properties). Quantum me- 
chanics does not predict the result of a single mea- 
surement, but guarantees (x) = lim 4r %i = 0. 

For each measurement with the outcome xt ^ one 
can say that P-invariance is broken in this particu- 
lar experiment, "event-by-event". In this simple case 
"breaking" is clearly of statistical origin and has noth- 
ing to do with dynamics - i.e. properties of the po- 
tential V(x). Therefore it is common in quantum me- 
chanics not to use such terminology and compute in- 
stead nonzero P-parity even observables, such as, e.g. 
(x 2 ) — J x 2 dx\ipo(x, t)\ 2 7^ 0, characterizing the pattern 
of quantum fluctuations. What is however important is 
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the textbook average over events / average over proba- 
bility density equivalence. 

By way of another simple analogy consider a system 
of massless quantum fields subject to boundary condi- 
tions at typical distance scale L characterized by a unit 
3-vector n. To be concrete one can think of electro- 
magnetic Casimir vacuum between parallel plates at dis- 
tance L with n being normal to the planes. Let this vec- 
tor smoothly fluctuates in random directions with typ- 
ical frequency w, which is assumed to be much smaller 
than c/L. One studies the quantum average of energy- 
momentum tensor for the fields, (T^ v (x)). Since the 
problem is quasi-stationary, the general answer is given 
by 

(Tfj,u(x)) = a(x) + b(x) n^n v + 0(ujL/c) (11) 

On the other hand, performing an average over time 
period T >• cj _1 one should have 



- J dt(T llv {x))=a{x)g llv (12) 



since no memory has remained about the particular di- 
rection the vector n is pointing to. Thus experiments 
with the detector time resolution ujT <C 1 will observe 
0(3) violating local answer (| 1 1 1) while those over long 
time scales ujT ^> 1 will see 0(3) respecting answer 
(|T2"j) . It is of crucial importance that some physical pro- 
cess with the typical life time scale comparable or larger 
than the plasma life time does exist and it is responsi- 
ble for creation of P- and CP-odd domains in the dense 
and hot matter in Minkowski space-time. It seems to 
be rather subtle point in this case how a relation be- 
tween Euclidean expression ([2]) and Minkowskian dy- 
namics should look like. In any case the existence of 
scale separation between the process dynamical scales 
and the measuring device ones is necessary for the whole 
picture to make sense. 

Since detailed picture of the discussed microscopic 
quantum/classical interplay is beyond us, our attitude 
here is purely phenomenological. We define the effective 
^-dependent current J^(x,rf) as 



function in r\ and 



(13) 



where electric current j^ix) = ^(x)Qj fJi 'ip(x) with 
quarks charge matrix Q = diag(2/3, — 1/3, — 1/3). The 
state \^n} is characterized by 



/ d 4 ydf(y)\n v ) =77 
Jv 



(14) 



dr]J^(x,r]) = 



(15) 



Since by assumption each event is characterized by some 
value of ?7, positive or negative with equal probabil- 
ity, this corresponds to "averaging to zero" over many 
events. 

To proceed further it is convenient to use the formal- 
ism of partial partition functions: 



z = J D^cM-sm\{ J dm ~8{m-Oi[$\ 



(16) 



where S[$] is the standard Euclidean QCD action, $ 
stays for dynamical quark and gluon fields A,tp,ij} and 
Oj[<&] is a gauge-invariant operator made of these fields. 
We approximate the real detector with finite resolution 
by the choice of the "detector function" S(x) in Gaussian 
form: 



^ = 27 



dAexp(-A 2 ;7 2 + ^) (17) 



so that dr]8(r]) = 1. 

We are interested in a value of the electric current 
(fT3|) . For exactly conserved axial current dj 5 = one 
would have (0|j M (a;) • dj 5 (y)\fl) = 0. Due to (electro- 
magnetic) anomaly however the result reads (for isovec- 
tor components) 

i [ dxe^-v) • dj 5 < a (y)\Q) = 



= Tr[QV]-^-^)-^ !/ (18) 

where t a are generators of flavour SU(2) or SU(3). 

For singlet current the anomaly gets gluon contribu- 
tion 



dj* = -Tr[Q 2 ] ^F^F» V 



N 



167T' 



Tr G^G^ (19) 



(notice that for uniform magnetic field F^F^" = 0) 
and computing 



where 

!iy = 



z 



D* j M (x) 5 (t} - n v ) exp(-S[$]) (20) 



d 4 ydf = - 



N 



f 



d 4 yTvG^G^ (21) 



at the leading order of the cluster expansion 
(A exp B) « (AB) exp((B 2 )/2) 



(22) 



It is physically obvious that J^x^) must be an odd 



valid for {A) — and (B) = 0, one gets in this approxi- 
mation 



J^r,) = -Tr[Q 2 }-^ 



N c V e-^/ 2L ' 2 



4tt 2 72^16 



d 4 q 
(2nY 



e lqx fv(q)iq" 



■ F 



(23) 



Here L 2 = I 2 + (n v ) and the formfactor is given by 
fv{o) = jy d V ex^>{—iqy). In the infinite volume limit 
X = lim (n v )/VN 2 where x is the standard topologi- 

cal susceptibility. 

The expression (12^1) deserves a few comments. First, 
the right hand side of (|23|) is odd function of r] as it 
should be, and at small 7/ the current is linear in r\. If 
the point x is far from y G V the current vanishes due 
to formfactor fv(q), i-e. the current flows only in the 
interaction volume V. On the other hand, if x 6 V and 
V is large enough to neglect surface terms, the current 
also vanishes as it should be for any finite- volume effect. 
The volume scaling (n v ) ~ V for the phase with finite 
correlation length is another manifestation of the same 
fact. 

It is worth mentioning that the maximal current is 
reached at r\ ~ L and decrease as J max cx BJ tL 2 (where 
r ~ V 1 / 4 ). This result seems counter-intuitive. Indeed, 
a naive picture would suggest that stronger fluctuations 
of topological charge (n v ) are to correspond to stronger 
currents J^{x,rf). This in fact is not the case. Rough 
physical explanation follows from (fT8|) : since the prod- 
uct of and dj 5 is fixed by electromagnetic anomaly 
(i.e. by the magnitude of external abelian field F^) 
large dj 5 corresponds to small j M and vice versa. Let 
us remind that according to the lattice data the 
magnitude of topological charge fluctuations experience 
rather sharp drop above the deconfinement transition. 
According to the above it means the effective enhance- 
ment of maximal possible electric current fluctuations! 
Of course at too small (fly) Gaussian approximation 
(neglect of higher order correlators) we have used is to 
break down. 

It is seen that the discussed effect is a result of subtle 
interplay between strong and electromagnetic anoma- 
lies (see related remarks in Q). While the later one is 
responsible for correlation between vector and axial cur- 
rents, the former anomaly provides non-conservation of 
axial charge due to topological nonperturbative gluon 
fluctuations. The question about /15 distribution ad- 
dressed in the introduction is translated here into the 
question about rj distribution for experimental events. 



IV. CHARGE FLUCTUATIONS ASYMMETRY 
AND POLARIZATION OPERATOR 

Perhaps the most logically consistent way is to study 
transition matrix elements of (JXJ) between states of op- 
posite P-parity. This corresponds to: 



where the states |f2) and \A) have opposite P-parities 
and (A\n) — 0. Of course the expression (|2"4")l is noth- 
ing but the electromagnetic polarization operator in the 
state |f2) saturated by particular states in spectral ex- 
pansion. 

This line of reasoning has been addressed in the liter- 
ature before. Local averages like (j 2 (x)) were computed 
in pioneering studies of CME on the lattice [H H|] and 
many interesting patterns were found. Later nonlocal 
averages {j l x{x)j v {y)) are computed [45|,[46j]. We find it 
worth reminding once again that since the typical cor- 
relators we are interested in are given by dimension six 
operators, their local matrix elements are strongly UV- 
singular 

(jl (i)) F «A 6 + F 2 A 2 + UV-finite (25) 

where A is UV-cutoff and F-external field strength. 
Even subtracted average (j 2 (x))F — (j 2 {x))o is diver- 
gent. This problem is overcome in numerical lattice 
calculations, but present analytical challenge for any 
attempt to describe CME in terms of local matrix ele- 
ments. To our view this is a clear signal about intrinsic 
nonlocal nature of the discussed phenomenon. 

Polarization operator in the CME context is studied 
in [22j . There are two main differences between our 
approach and that of the cited paper. First the regu- 
lar contribution (given by polarization operator in mag- 
netic field) and CME-contribution ( pro portional to ^5) 
are separated from the beginning in [22j (in some sense, 
quantum currents are superimposed on top of the clas- 
sical current (1)). Wc follow another logic and consider 
polarization operator as the only source of asymmetric 
charge fluctuations, but extract a particular formfactor 
from it, which corresponds to negative parity interme- 
diate states. Second, the expression for charge fluctua- 
tions observable as a functional depending on polariza- 
tion operator is different in our paper from that of [22| . 
We will make more comments on that below. 

In this section we discuss P-odd x P-odd contribu- 
tions to P-even observable, the role of which is played 
by current correlator (jfj,j v )- It seems physically clear 
that this object should contain some information about 
charge distribution The exact form of this corre- 
spondence is however far from trivial. One could think 
of several ways to relate these quantities. Before pre- 
senting our approach to this problem let us mention 
other methods used in the literature. First, we notice 
that the current in ^-direction is given by 



(n\jij k \n) -> Y,m\ A )( A \jk\n) (24) 

A 



e z j3 + e y j2 =yil+ if i e z sin + e y cos 4>) (26) 



G 



and the corresponding charge difference from ([JJ is 

«*(*+-*-),, r (27) 



where by the brackets (...) e we denote the average over 
events. One has ((a + — a_) 2 ) e oc (j'f + jf ) where the 
current product is assumed to be local. This is very 
close (but different) to the definition used in 43]. It is 
natural to expect that positive definite ((a+ — a_) 2 ) e 
should be nonzero even without an y m agnetic field. 

Another relation is suggested in [22|]. It is written in 
terms of event average of the cosine, where a, /3 = +, — 
and N± is the total number of outgoing particles of a 
given charge: 



(COS(0 Q + <j)p'))e OC 



a/3 
N a N p 



(J 



■ft) 



(28) 



where, up to some background terms 

(cos(^ Q + <j>p')) e = (vi, a Vi,p) e - {a a af}) e (29) 

Assuming charge independence of vi <a and equal num- 
bers of particle species N + = N- = N one gets 
((o+ — a-) 2 ) e oc Of — j|) if one neglects v\ >ot term with 
respect to a a term. In fact, the leading term, which is 
always contained in jg component, coincides for both 
expressions, while the procedure of taking into account 
fluctuations in the reaction plane is different. 

In this paper we use alternative signature provided by 
charge density fluctuations and not spatial components 
of the currents. An attractive feature of this quantity 
is that it is well defined even in the static limit. To this 
end consider electric charge in some spatial volume V 
at temperature T: 



eQv = e / dxj (x) 
Jv 



(30) 



Since we work in zero density approximation the quan- 
tum average of this object vanish: 

(Qv) = (31) 

This is not the case for its square: 

(Q^)=-k[ dx [ dx' 1144(2, x') (32) 
Jv Jv 

where 1144(0;, x') is Euclidean polarization operator in 
constant external field F^ v and at temperature T 
Wick-rotated from the standard Minkowski expression 

II 00 (x,x')\ 

IlW (*,*') = i{T{3 li {x)j v ^)})F,T (33) 

with = 'tpQ^fiij); Il^f' ■<-» II^P, notice the sign con- 
vention (|3"2")l corresponding to positive-definite (Qy) in 
the static limit. In the standard way we denote 

IV (9) = / d^x e -*K»-*') n^(x, x') (34) 



with fj,, v = 1,2,3,4 and q = (91,92,93), q± = (91,92)- 

The operator k in (j5"2l accounts for temporal profile of 
the process. In terms of momentum space components, 
(l3"2"j) takes the following form 



(Qv) 



dfU 
2tt 



K(q 4 



dq 



|^(q)| 2 n 44 (q,94) 



(35) 

where the form- factor i*V(q) = J v cixexp(iqx) keeps in- 
formation about spatial profile of the volume V, while 
the temporal factor ^(94) = / dTg(r) exp(zg4r) encodes 
temporal (in Euclidean sense) profile. For finite temper- 
ature case we consider here the standard Matsubara re- 
placements 94 oj n = 2i:nT and (27r) _1 J dq 4 — > T^2 n 
are to be performed. The choice g(r) — T we will adopt 
in the rest of the paper physically corresponds to the 
static limit where only the lowest Matsubara frequency 
n = contributes: 



(Q 



VI st 



= -T 



dq 



|Fy(q)| 2 n 4 4(q,0) 



(36) 



Using the expressions from Appendix it can be checked 
that in thermodynamic limit V — > 00 without external 
field one reproduces standard Stcfan-Boltzmann answer 
for elementary fermions lim (e 2 Q v ) s t/V = e 2 T 3 /3. In 

V— >oo 

case of quarks one should of course understand eB as 
qjeB and introduce additional trace over flavors with 

the factor N C Q 2 : ^f A ' T -> A^ c J2f q)^®^ ■ For the 
sake of brevity we will use the simple notation as for 
elementary fermions of unit electric charge having in 
mind the necessity to make the replacement discussed 
above in the final answers. 

In the limiting case of no background B — 0,T = 
one has IL^q, 94) = q 2 II(9 2 ) and, at the leading order, 
for large 4- volumes V4: 



(<2y>B=o,T= oc n'(o) • v 4 ~ 1/2 



(37) 



where the condition of gauge invariance 11(0) = has 
been taken into account and the volume V4 = R 3 x t 
is assumed to be uniform: R ~ t. Thus the expression 
(|32|) is UV-safe and vacuum charge fluctuations in a 
given space-time region is purely finite- volume effect. 

We can now come back to the definition (|32]l and 
rewrite the coordinate integration in cylinder coordi- 
nates with the axis 1 as the polar axis and angle <p de- 
fined in the 23 plane. This is the same notation as in 
notice that in the standard setup azimuthal angle is 
usually defined in the plane 12. This allows to represent 
the form- factor Fv{q) as 



dx\ e 



Rdp 
Jo 



(38) 



where qp = 92X2 + 93X3 = 92/9 cost/) + q^p sin</> and the 
structure of integration upper limit is determined by the 
chosen model for spatial distribution (sharp boundary, 
smoothed boundary, Gaussian shape, exponential shape 
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etc). The sin</> - mode in Fourier expansion of (|38)) is the change -FV(q) — > fv {qi, ?2, 93) where 
multiplied by the following coefficient 



ci = (1/tt) 



> sin ©e 



2iq 3 



' Y J|il,,i (39) 



fv(qi,Q2,q3) 



2iq 3 



dx x e lqiXl / pJi{qp)dp (41) 



, - AIT^ 2 Ti, i r r In the same way Uq v ) ) corresponds to the exchange 

where q — \ qi + qi. lhus we have for expansion of , . , v , „ , . . m trr**. , 



in harmonics: 



sin^ / #'sin0'((q^) 2 ) + 



(0? 

./ n ./ n 

( 4 °) 

where ((q v ) 2 ) is given by the same expression (f36| with 



93 <Z2 and sin <p f-> cos 0. Making use of ((4]) , (|27l) and 
(|40[) we obtain the following relation for the asymmetry 

(q 2 v) = ((q v ) 2 ) - ((qv 1 ) 2 ) = £ «/3cos(0 a + V) 

(42) 



(4) = iV 2 -(((a + -a_) 2 ) e -<K 



) 2 >e)=^ 



f/q 



(2tt)3 gf 



qi 



dxx e lqiXl / pJi(qp)dp 



n 44 (q,0) (43) 



It is obvious that the above expression has to be pro- 
portional to magnetic field since there is no other 0(3)- 
violating factors in the problem. The effect we are look- 
ing for corresponds to strong enhancement of (|43[) in 
external magnetic field and hence, from experimental 
point of view, strong dependence of (|43[) on centrality. 
It is to be stressed that the multiplicity factor N 2 is by 
itself strongly centrality-dependent. This dependence is 
kinematical and has nothing to do with magnetic field 
dependence of (q v ). Only the latter lies at the heart of 
CME. 



V. GENERAL STRUCTURE OF 
POLARIZATION OPERATOR 

In this section we analyze general structure of polar- 
ization operator in the background of nonzero temper- 
ature and magnetic field. As is clear from the above 
discussion, this is a necessary prerequisite before one 
can compute the charge fluctuation asymmetry (14*31 . 

First of all let us make a few general comments about 
space-time dependence of current-current correlator. In 
confinement phase (i.e. at sufficiently low temperatures) 
at large distances and for weak magnetic field one ex- 
pects general structure of Euclidean polarization oper- 
ator of the following form 



(j{x)j(x')) oc e-™^-*'! + C{B) ■ e ~ m ^ x ~ 



(44) 



with C(B) cx B 2 . This interesting effect of different 
parity states mixing in external field is similar to the 
one observed long time ago in 47] at finite tempera- 
ture. The long-distance correlations are thus saturated 
by the lightest degrees of freedom (i.e. pions in the con- 
finement phase). On the other hand, in deconfinement 
phase at strong fields, if Larmor radius is much smaller 



than Aqp^, no quarks can propagate in transverse di- 
rection at all: 



(j(x)j(x')) oc 



-eB(x-x') 2 ± /2 



(45) 



Large- iV c suppressed transverse correlations are possible 
only due to gluon degrees of freedom. 

We confine our attention in what follows to a particu- 
lar case of purely magnetic constant abelian background 
field in the thermal bath rest frame at nonzero tem- 
perature T. We have chosen F i2 = — F21 = B, i.e. 
magnetic field is directed along the third axis. The 
temperature effects break Lorentz-invariance and the 
physical answers depend on 4-vector which repre- 
sents four- velocity of the thermal bath. It is normalized 
as u^u^ = 1. In the present paper we take zero chemical 
potential p — 0. It is to be noticed that many general 
conclusions concerning the structure of polarization op- 
erator stay intact for p 7^ since the latter is associated 
with the same four-vector given by — (1,0,0,0) 
in the medium rest frame. 

The polarization operator (|34[) is a rank two tensor 
depending on two polar vectors g M and and anti- 
symmetric tensor F^ u . The general decomposition of 
(I34[) in terms of independent tensors was extensi vely 
studied in the literature starting from [HI, H^] , see [5C| 
for recent exposition and [51( for a useful collection of 
references. Generally, one is to deal with 4x4 = 16 
independent tensor structures, built by multiplying the 
four independent base vectors g^, u^, q a F afl , q a F@ Fp^. 
It can be shown however that general requirements of 
being transversal 



q^ u (q) = q v U^{q) = 



(46) 



and Bose symmetric Il^^g) = 11,,^ (— q) together with 
generalized Furry's theorem [48j 



ILj^g, u, F) = II^(g, -u, -F) 



(47) 
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reduce the number of independent tensor structures to 
six. Two of them are field- independent, the other two 
depend on linearly and the last two - quadratically 
(notice that our numeration of the tensors is different 
from the one adopted in [48| ) . Their explicit form reads 



= (<7 2 % - 9m ( u <7)) (<l 2 u v - q v (uq)) 



$(3) 

(IV 



$(4) 

(IV 



(uq)(q p F vp q p 



q v F w q p + q 2 F llv ) 



(up,F vp q p 



u v F w q p 



(uq)F pv ) 



lj>( 5 ) — P nPP ri 



*$ - (<? 2 <W - q„Q P )FSF aa (<fS av - q a q u ) (48) 
The coefficient functions of the decomposition 



IL^, u (q,u,F) 



(49) 



depend on gr 2 , mixed invariants (uq) 2 , (qF) 2 , (uF) 2 , 
(qFu) 2 , pure field invariants F 2 ,FF and also the tem- 
perature T and particle data, encoded in matrices Q 
and M. The expression (|49|) allows to discuss current 
correlations asymmetries in invariant way in any theory 
where the expression for polarization operator can be 
obtained. 

Having these general prerequisites let us come back to 
analysis of correlation patterns. For our choice F12 = B 
the invariants (uF) 2 , (qFu) 2 and FF equal to zero. In 
what follows we will be especially interested in a partic- 
ular type of contribution to Ii pl ,(q) proportional to the 
(7) 

tensor structure given by the product of two axial 
vectors 



= p 



(50) 



It is not independent and one easily checks that ^ p 
can be expressed as a linear combination of 



(7) 



qH$ = (q 2 F 2 /2 - (qF) 2 ) *£> 



MS (51) 



Let us consider tensor structure of the polarization 
operator in more details. First of all, since we are inter- 
ested only in diagonal 11, 22, 33, 44 and also 34 compo- 
nents in this paper, we have no contributions from 



and tt^ 4 ) because the tensors and \&}J are antisym- 
metric and also vanish for fi = 3, v = 4 in the chosen 
background field. Second, we notice that for /1, v equal 



to 3 or 4, one has identically SS>$ = 0. Adopting con- 
ventional notation: q± = ((71,^2), Q\\ — (93,94) we can 
rewrite ((49]) using (|51"]) as 



(1) 



(52) 



II||(g) = 7T ■ *|| -r „ ■ *|| 

where the new invariant functions are given by 

7T»> = 7T« - (q 2 F 2 12 - (qF) 2 ) 7r< 6 > 



n (T) = ^(2) . -(F) = q 2 w (6) 



(53) 



As for the diagonal correlators in 12-plane, one has 

U ± (q) = ttW) • + 7r< T ) ■ *f + ^ (F) ■ *f (54) 

where and ir^ are defined by the same expressions 
(l53l) while tt^ form-factor reads 



(55) 



It is seen that the correlators of our interest can be 
decomposed into just three independent structures. The 
first, corresponds to purely quantum fluctuations. 
It has nonzero limit at both B — > and T — > 0, which 
coincides in this case with the textbook expression for 
polarization operator. The second structure, ir^ is re- 
sponsible for thermal fluctuations. It vanishes at T — > 0. 
It is worth mentioning that both functions n^' and 
7r( T ) depend on temperature and external field (since 
the pattern of both quantum and thermal fluctuations 
is sensitive to the external conditions) and our notation 
corresponds rather to the limiting form of these func- 
tions. 

We notice that the terms proportional to and 
7r( T ) are identical in (j52"j) and ([53]) up to obvious change 
of notation || ^5 _L. This is to be expected since quan- 
tum and thermal fluctuation are 0(3)-isotropic. The 
only non-isotropic term (and the most interesting for us 
here) is the last terms: in (|52|) and in ([54]) . 

The former one takes into account charge (and also the 
current component j'3) fluctuations induced by the ex- 
ternal magnetic field. P-parity structure of this term is 
given by 



P-even = P-evcn 



F 3p f 
axial 



axial 



It is to be compared with the thermal contribution pro- 

(2) 

portional to : 

M.73.73) = 7r(T) x P3(up) x p 3 (up) 
P-even = P-even x vector x vector 

This directly corresponds to our discussion in the in- 
troduction: in the latter case the thermal fluctuations 
are distributed isotropically in the thermal bath rest 
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frame, while in the former one there are electric currents 
fluctuating along the magnetic field. The magnitude of 
these fluctuations is measured by the function tt^ f > , and 
no physical principle forces it to vanish either below or 
above critical temperature. Physically corresponds 
to P-odd intermediate states in the polarization opera- 
tor. 

The function tt( f > entering (1541) is a sum of two terms 
according to (|55|) . This also is to be expected. Charged 
particles flowing in the plane perpendicular to the mag- 
netic field are deflected by the Lorentz force, and this 
diamagnetic effect is taken into account by the form- 
factor 7r( 5 ) . It is absent in ILi. But particle's spin in- 
teracts with the field by means of cr Q ^F Q ^ term in 
as well as in IT_l which results in the factor gV 6 > in 
both expressions ([52]) and (f54|) . It is worth noting that 
according to our general logic the electric charge asym- 
metry is computed for the full expression for II44, not 
just from some part of it, proportional to . Thus it 
is legitimate to speak about CME-interpretation of the 
answer (|43[) only in the limiting case when provides 
dominant contribution. We discuss that in more details 
below. 



VI. MODEL EXAMPLES 

We analyze in this section two limiting cases where 
one can construct in explicit way. The first one 
corresponds to weak magnetic fields in the confinement 
phase. In this case the intermediate states are hadron 
resonances of negative P-parity (see closely related dis- 
cussion in [52j). However to select explicitely physical 
states making dominant contribution is far from triv- 
ial and the answer strongly depends on kinematics. We 
confine ourselves in this paper to the simplest case keep- 
ing only three neutral ^ intermediate states: 7r°, 77, 77'. 
Technically it is more convenient to consider from the 
very beginning matrix elements of vector currents be- 
tween vacuum and these states in external field. Mak- 
ing use of the definition of off-shell vector-vector-axial 
form-factor = J r 7r o* 7 , 7 , (q 2 , q 2 , q 2 ) (see, e.g. [53j]) 
with q = qi + q 2 



dx I dye 



iqix+iq 2 y I 



\Tr{Mx)j v (y)}\n°(q)) 



= e tlvafj q^q^J : 7! {q 2 , q\, gf) (56) 

one gets at the leading order in constant external field: 

(0\j^-q)W°(q))F=ieq p F w ^(q 2 ,q 2 ,0). (57) 

The expressions for rj and 7/ contributions are com- 
pletely analogous with the replacement of T v by and 

Thus the g 2 -dependence of polarization operator in 
external field is determined in this approximation by 



the form- factors J r ( / > (q 2 ,q' 2 , 0) with one on-shell leg (cor- 
responding to external field vertex) . These form- factors 
are essentially nonperturbative QCD objects. Let us re- 
mind that on-shell (i.e. at the point ^(m^, 0, 0)) they 
are fixed by triangle anomaly, for example for pion: 

^(^,0,0) = -^^- (58) 

Another important case is large q 2 — > 00 limit where one 
has (for chiral fermions) J-^iq 2 , q 2 , 0) — !• XfF w /3 where 
Xf is QCD quark condensate magnetic susceptibility, 
defined by (0\qa ^q\0) f = e q XF{qq)F f _ LU . Different ap- 
proximation schemes valid at intermediate momenta are 
discussed in the literature (see, e.g. [Hj]). 

Having written the field-dependent matrix element 
([57)1 one is able to express the invariant function 



(q ) as follows: 



l^(<zW,Q)l ; 

q - ™>i 



(59) 



From the point of view of expression (|43[) the domi- 
nant contribution to asymmetry is this phase comes 
from the lightest degree of freedom, i.e. massless in 
the chiral limit pion (to be more precise, we assume 
the limit m^R -C 1). Choosing for concretness Gaus- 
sian boundary condition (i.e. introducing the factors 
exp(— <jf i? 2 /2) into (|4"3"j) one obtains 



(q 2 v 



>=7 



6 -^) Ti? 3 

Fir 



(60) 



where the numerical factor 7 = 1.6-10 is of course spe- 
cific for this boundary choice. Certainly the result triv- 
ially follows from dimensional considerations. We see 
(qy) <§; 1 for phenomenologically reasonable choice of 
parameters. Contributions of mass gapped states bring 
additional suppression (and, in particular, break ~ R 3 
scaling). 

As the second example we consider free fermions in 
strong field limit. This regime would correspond to de- 
confinement phase where proper dynamical degrees of 
freedom are quarks and gluons with perturbatively weak 
interaction between each other. To compute polariza- 
tion operator under external conditions in perturbation 
theory one usually makes use of Schwinger proper-time 
technique and there is extensive literature on the subject 
[5"5U59l | where different kinds of external backgrounds 
were studied. The polarization operator in constant 
magnetic field and at nonzero temperature was calcu- 
lated in (6pj in imaginary time formalism. Our aim here 
is to put these results in a charge fluctuations asymme- 
try prospective. For the reader's convenience we repro- 
duce the explicit one-loop expressions for polarization 
operator IT| | given by [60( in Appendix of the present 
paper. 

It is convenient to present the Euclidean polarization 
operator in the following form 



IW<?±,<73,™) 



J2 A M^ Hq) + QM (61) 
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where the sum includes integration over proper-times 
and summation over Matsubara frequencies (see expres- 
sion (|A.3|) in the Appendix), the functions -A^g] poly- 
nomially depend on momenta components q, and the 
universal Euclidean phase <p(q) is given by expression 
(|A.4|) . The contact terms Q^ v {q) have no sensitivity to 
infrared parameters (like temperature or external field) 
and provide correct limit of II M „ at vanishing back- 
ground. 

One can notice that tv- f ^ can be simply related to 
the polarization operator components. Namely, solving 
the system of three linear equations (f52"j) for the choices 
(/if) = 44, 33 and 34 one finds all three invariant form- 
factors, including 



2nTn. 



B 2~(f) _ g3g4n 44 + (gj + gi)n 34 



(62) 



where q\ = q\ + q\ and q\ = u 

Thus the chiral magnetic form-factor is a nontrivial 
linear combination of II34 and II44. First of all we are 
to check that at B — > the r.h.s. of (|62|) vanishes. This 
is obvious at zero temperature since in this case there 
is the only tensor structure given by and 



93<74*44 + ill + 93)*34 = 



2WrM - , 



(63) 



For temperature-dependent parts it is rather nontrivial, 
the proof that this is indeed the case can be found in 
Appendix. 

The explicit expression for 7f( F ) looks especially sim- 
ple in small T regime. It reads 



7^) = 



1 1 

(4tt) 2 



du I dv ((l-w 2 )cotht2 + / ± (u,w))exp(-(/) (0) ) 



(64) 



where u — ueB and the functions <jr- > and fj_(u,v) are 
given in the Appendix. Notice that such form-factor 
was discussed in a different context in (49| . 
In the weak field limit one has 



hm ^ 

B->0 



1 

6^2 



1 dv C-' 2 )P-" 8 ) a (65) 
1 (4m 2 + (1 - v 2 )q 2 ) 2 



In the strong field limit (still at small T) the situation 
becomes more interesting - form-factor 7r( F ) provides 
dominant contribution to the polarization operator: 



n 



44 



qiieBf^ 



eB 



dv 



(1 - v 2 )q 2 



Am 2 + (1 - v 2 )ql 



(66) 



up to the terms O (g^/ei?) . One can say that all asym- 
metry of charge fluctuations is due to CME-like form- 
factor in this limit. 

We see another interesting effect - in the chiral limit 
does not depend on g 3 at all, while the dependence 
on q± is suppressed by the field B. On the other hand, 
the essence of the asymmetry of interest is just differ- 
ent dependence of the polarization operator on different 
components of momentum. Since the polarization op- 
erator itself linearly rise with B for strong field it is nor 
a priori clear which effect is to win. Detailed calcula- 
tion shows that in fact they balance each other and the 
asymmetry (1431) is not asymptotically rising with B - 
there is an effect of saturation. It is reasonable to sepa- 
rate different regimes depending on ratios between basic 
parameters such as B, m, T and R where the latter one 
stays for the typical 3-dimensional size of the volume V3. 



For two light flavors one can safely neglect quark masses 
m. Three other parameters are in the ballpark of 100 
MeV (for large fireball one can think of phenomenolog- 
ically realistic eBR 2 =5-7- 10). Without intention to 
cook up numerical factors but just to get feeling of the 
numbers, pluging (l66l) into (|43|) we get 



(q 2 v ) = 7' • RT 



(67) 



where again the numerical factor 7' = 4.1 ■ 10~ 2 corre- 
sponds to Gaussian boundary shape. Thus for asymp- 
totically large B one reaches " kinematical limit" for the 
asymmetry in our picture, despite numerically it is still 
very small. 



VII. CONCLUSIONS 

We have discussed three possible ways to study quan- 
tum physics behind chiral magnetic effect and elec- 
tric charge fluctuation asymmetry observed in heavy 
ion collisions. For all approaches the importance of 
scale separation is stressed - there should by hierar- 
chy of dynamical scales characterizing the life of quark- 
gluon phase after the collision and intrinsic QCD scales 
(perhaps field/temperature shifted) characterizing the 
nonabelian topological charge fluctuation pattern. The 
physical essence of CME as we tried to present it here is 
that the quark-gluon medium plays the role of a measur- 
ing device with respect to the topological QCD vacuum 
with the final particles electric charge asymmetry as an 
outcome. This is most clearly illustrated by the expres- 
sion (|2"3"j). 

The third approach we have considered, i.e. the anal- 
ysis of P-odd x P-odd contributions to P-even observ- 
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ables, is somewhat different because it provides nonzero 
results even for free fermions in magnetic field, i.e. with- 
out any "topological origin". We believe that this can 
be considered as a particular case of CME as well. Just 
nonzero matrix element of the vector current between 
vacuum and J ^ states in external magnetic field leads 
to asymmetric charge / current pattern as if there is fluc- 
tuating vector current collinear to B. Of course the de- 
tailed picture depends on the actual quantum dynamics 
of these J h degrees of freedom, and we have shown 
that indeed it is strongly suppressed in the confinement 
phase. Nevertheless we find it legitimate to interpret 
this dynamics using the same CME-like language since 
namely this anomaly-driven vector-axial correlation is 
at the heart of the effect, while the concrete way of life 
of the axial degrees of freedom (distribution function 
for /15 in the standard CME analysis) is of secondary 
importance. 

We have left without attention all aspects of temper- 
ature dynamics in this paper. Despite no drastic quali- 

I 



tative effects are expected it is interesting to study the 
asymmetry in the whole parameter space spanned by 
(B, T, m, R). This could clearly have phenomenological 
applications to heavy ion collision physics whose under- 
standing is the main challenge for modern QCD. 



ACKNOWLEDGMENTS 

The authors acknowledge useful discussions with 
P.Buividovich, A.Gorsky, D.Harzeev, A.Kaidalov, 
M.Polikarpov, Yu.Simonov and V.Zakharov. The work 
of one of the authors (V.O.) is supported by the "Dy- 
nasty" foundation personal grant. 



Appendix 

For the reader's convenience we present explicit ex- 
pressions for one-loop polarization operator as com- 
puted in [6(| ■ It reads 



^ n ' l=-oo 



(A.l) 



where Si(p) is fermion propagator in external constant 
magnetic field and Q^ v (q) is the " contact term" needed 
to cancel ultraviolet divergencies. It has no depen- 
dence on soft backgrounds like temperature or exter- 
nal field. The sum goes over fermionic Matsubara fre- 
quencies lji = (21 + 1)ttT, while the bosonic one is 
q± = ui n = 2nirT. 

Thus the general structure of Euclidean polarization 
operator is given by 



IV(4_L,*3,n) = ^A^(q) 



+ QM (A.2) 



The sum is given by the following expression: 

V Je J ~ 1 l=-oo 

The universal phase 4>{q) has the form 

m = (o) (<z) + uwf 



q 2 ^ cosh u — cosh uv 



eB 



2 sinh i 



m 2 + W, 



(A.4) 



where Wi = Cji — ^^-ui n and u — ueB. The contact 
term is given by 

00 

„ 1 f du 



12tt 2 



(<? 2 <W - q»q v ) (A.5) 



r 



The function A^ v (q) polynomially depends on momenta 
components q and for 3,4 components reads 



^44(9) = cothw 



l 



2Wf + q 4 vWi 



1 



"?3 



A^{q) = q 3 vWi + 
^33(9) = -cothit 



1 



-Q-i 



<&- 



cothu 



q^vWi 



(A.6) 
(A.7) 



-f±(u,v) 



where 



fx(u,v) 



v coth u sinh ? 



cosh? 



sinh it 



(A. 



(A.9) 



To get zero temperature limit of the above expressions 
Poisson summation formula is useful 



£ < 

/— — 00 



— 2-nizk 



(A.10) 



k= — 00 
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In particular one gets 



oo 1 

lim T V e- uW ' = —= (A.ll) 

l=-oo V 



Another necessary demonstration of self-consistency 
is a proof of vanishing of B 2 tt^ defined by (j()2"j) at 



B = for any T. One has, at B -> 
9394^44 + q 2 v4 34 = ^ U Q - + uW^ 2 ) = 

= » d ^ -^cwf+^A (A . 12) 

It is easy to check that the latter expression gives zero 
result when integrated from -1 to 1. 
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